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Coupling of transverse and longitudinal response in stiff polymers 



O 

o 

On: 



o 

c3 



i 

c 

o 
o 



(N 
> 

m 

(N 

m 

o 
o 



X 



Benedikt Obermayer^B and Oskar Hallatschek 2 '[l] 

1 Arnold Sommerfeld Center and Center for NanoScience, 
Ludwig-Maximilians-Universitat Miinchen, Theresienstr. 37, 80333 Miinchen, Germany 
2 Lyman Laboratory of Physics, Harvard University, Cambridge, MA 02138, USA 

(Dated: February 1, 2008) 

The time-dependent transverse response of stiff inextensible polymers is well understood on the 
linear level, where transverse and longitudinal displacements evolve independently. We show that 
for times beyond a characteristic time tt, longitudinal friction considerably weakens the response 
compared to the widely used linear response predictions. The corresponding feedback mechanism is 
explained by scaling arguments and quantified by a systematic theory. Our scaling laws and exact 
solutions for the transverse response apply to cytoskeletal filaments as well as DNA under tension. 
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In tracing back the viscoelasticity of the cell to prop- 
erties of its constituents, a detailed understanding of 
the mechanical response of single cytoskeletal filaments 
is indispensable. Due to their large bending stiffness, 
these filaments exhibit highly anisotropic static [l[ and 
dynamic 0, H, 3 features, such as the anomalous in- 
growth of fluctuation amplitudes in the transverse direc- 
tion 0,0, i.e., perpendicular to the local tangent. The 
related response to a localized transverse driving force 
has so far been examined only by neglecting longitudi- 
nal degrees of freedom Q, although these polymers 
are virtually inextensible, and transverse and longitudi- 
nal contour deformations therefore coupled. In this Let- 
ter we show that longitudinal motion strongly affects the 
transverse response even for weakly-bending filaments 
and leads to relevant nonlinearities beyond a character- 
istic time if. 

The physical key factors controlling the transverse re- 
sponse may be understood from Fig. [TJ which shows a 
weakly-bending polymer (bending undulations are exag- 
gerated for visualization) shortly after a transverse driv- 
ing force f± has been applied in the bulk. In response 
to this force, the contour develops a bulge. Due to the 
backbone incxtensibility, this bulge can continue grow- 
ing only by pulling in contour length from the filament's 
tails. This effectively reduces the thermal roughness of 
the contour d, 0, [l(| , at a rate substantially limited by 
longitudinal solvent friction. The resulting coupling to 
the longitudinal response tends to slow down the bulge 
growth. In order to describe this feedback mechanism, 
we start with a scaling analysis and treat the simpler 
athermal case first. To connect to the biologically im- 
portant situations ofprestressed actin networks and 
prestretched DNA II 211 , we then extend a recent theory of 
tension dynamics [13[ to calculate the nonlinear response 
for unstretched and prestretched initial conditions. 

Consider the overdamped dynamics of an initially 
straight stiff rod of total length L. Suddenly applying 
a transverse pulling force f±, for simplicity in the cen- 
ter of the rod, leads to the growth of a bulge deforma- 



tion. The generated friction in the transverse and lon- 
gitudinal direction needs to be balanced by correspond- 
ing driving forces. Viscous solvent friction is modeled 
via anisotropic friction coefficients (per length) £j_ and 
0i = ff_L with £ sa \ for transverse and longitudi- 
nal motion, respectively After a time t, the resulting 
bulge has some characteristic height A±(t) and width 
l±(t). The transverse force f± balances the drag force 
C±(± Aj_/i acting on a polymer section of length £± mov- 
ing transversely with velocity A±_/t through the solvent; 
hence, Aj_ ~ f±t/(£±£±). Naturally, the contour length 
along the deformed rod section is larger than its longitu- 
dinal extent £±. Assuming a simple "triangle" geometry 
as in the blow-up in Fig.[TJ the difference is roughly given 
by Aj_/^j_. In order to provide this stored (or excess) 
length, the filament's tails are pulled in by a longitudi- 
nal force /y. The latter has to balance the longitudinal 
friction that acts on the filament's tails of length L mov- 
ing longitudinally with a velocity given by the tempo- 




FIG. 1: (Color online) A transverse point force f± applied to 
the contour r(s) (dark) translates, through the formation of a 
bulge of height Ax and width i±, into a longitudinal pulling 
force /|| acting on the polymer's tails. This force induces back- 
bone tension f(s) (light) that penetrates the contour within a 
region of size £« where thermal undulations are straightened. 
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TABLE I: Summary of crossover scaling laws for an initially unstretched filament. The crossover time tf is implicitly de- 
fined through tf = /|j~ 2 (tf), /||(t) is the induced longitudinal force, Ax(t) is the transverse response and £±/\\ is the trans- 
verse/longitudinal correlation length @, H, O- ( a ) Athermal case, fn pulls in the filament's tails of length L. tf — (7o/±)~ 2 
with 70 = (Ci) 2/3 /I /3 - (b) Thermal case. The filament's tails have effective length t\\ -C L. U = (7/x) 2 with 7 = {(,£ P ) 2 ^' f±' ■ 



ral change of the excess contour length contained in the 
bulge. Estimating /|| ~ (\\LA 2 L /{£±t), we plug in Ax 
from above and get /y ~ (\\Lf1t/{(j_£ 3 ± ). 

The yet unknown time-dependent width £±{t) of the 
bulge is controlled by the relaxation spectrum of bend- 
ing deformations. In the weakly-bending limit, the trans- 
verse displacement field rj_{s,t) of an overdamped inex- 
tensible rod with bending stiffness k obeys [B| 

Cx$rL = -«rr + /,|(tK, (1) 

in the presence of a longitudinal pulling force f\\{t). 
Primes denote derivatives with respect to the arclength 
coordinate s £ [— 4, k\. In the following, we set n and 
Cx to unity, such that time is a length 4 and force a 
length -2 . From a simple scaling analysis of Eq. (JTJ) , 
r±/t ~ r±(£j 4 + f\\£± 2 ), we deduce the growing size 
£x(t) of a bending deformation (assuming £± «C L). In- 
serting appropriate formulas [l3( for £±{t) into the re- 
lations for Ax and h derived before finally yields the 
selfconsistent scaling laws for /y (t) and the nonlinear re- 
sponse Ax(t) summarized in Table QJa). For short times 
the coupling effect is irrelevant and Ax(t) is linear in 
/x- However, this requires the small force /y to pull 
in more and more contour length from the tails and in- 
creases the longitudinal friction to be balanced by fn . At 
the crossover time tf, this force becomes large enough 
(typically, f\\ ~ 7/x > /x) to feed back onto the trans- 
verse dynamics, which is manifest in nonlinear dependen- 
cies [18[ on /x. In particular, it considerably slows down 
the bulge growth, which in turn requires fn to pull in 
contour length at a slower rate and eventually makes it 
decrease. 

The essential difference for nonzero temperatures is 
the presence of thermal contour undulations, see Fig. [IJ 
which are correlated over the persistence length £ p = 
(fceT) -1 , and straightened out by the longitudinal force 
/|| . Still counteracted by longitudinal friction, this hap- 
pens first only within a small but growing region of size 
£\\{t) (see Refs. d !, ©, M, EDI)- Correspondingly, the 
force f\\(t) from above has to be generalized to a ten- 
sion field f{s, t), which decays over the length scale £\\ it). 
Crossover scaling laws for £\\{t) 1 shown in Table |jb), 
were derived for constant external force in Ref. [IH and 
can be generalized to (weakly) time-dependent "exter- 
nal" forces such as /||(t). The thermal problem is es- 
sentially analogous to the athermal case for late times 



t > t L where t L is defined via £\\{t L ) = L. However, if 
the region £\\{t), where the contour straightens, does not 

yet extend to the filament's ends {£\\ <C L, or t <C 
the "thermal" rod has only an effective time-dependent 
length of £\\{t). Hence, scaling laws for the nonlinear re- 
sponse are then obtained simply by replacing L — > £11 
111 Table Ufa), which gives the results summarized in Ta- 
ble Hljb). These apply to initially unstretched filaments 
while the general case of prestretched initial conditions 
is discussed below and summarized in Fig. [3] Naturally, 
the replacement L — > in affects only the long-time scal- 
ing of the nonlinear response Ax(t) - on short times 
f C (f, the transverse dynamics evolves undisturbed 
by the longitudinal one. We expect the anomalously 
slow long-time response to be observable in many bio- 
logical situations. In aqueous solution, we roughly es- 
timate a crossover time tf ~ I0~ 2 s/ fjJpN] 8 / 3 for typ- 
ical microtubules with I » 10 /^m [14| (representing 
the athermal case). Under thermal conditions, where 
the "interesting" time window is between tf and t^, we 
get t f ks lCT 3 s/7x[pN] 16 / 7 and t\ w 0.2s//x[pN] for 
(unstretched) actin filaments of about 20 /im length 0] , 
which implies that the actin response to myosin mo- 
tors becomes nonlinear on time scales comparable to 
the duration of a single power stroke [lB] ]. Filaments 
in actin networks (mesh size £ ~ j^L « 0.5 /im) un- 
der stresses of about 1 Pa arc usually so short that 
tf 3> 1 1 ~ 10 -4 s, but the coupling nonlinearity should be 
observable in the viscoelastic response Finally, tf « 
10~ 5 s//x[ P N] 16 A and 4 « 0.05 s/(/x[pN]/ prc [ P N] 5 /s) 
for DNA (L « 20 /jm [lj) prestretched with j prc < f±. 

In order to support and quantify the scaling picture 
developed above, we proceed with a systematic approach 
similar to Ref. [13[ based on the length scale separation 
£\\{t) ^> £±it). As long as the dynamics induced by the 
transverse force is not influenced by end effects {in <C L), 
we consider a semi-infinite arclength interval, s G [0, 00), 
and represent the transverse force as a boundary con- 
dition at s = 0. In the wormlike chain Hamiltonian, 
H = I Jds[r" 2 + fr' 2 ], the tension f{s,t) enforces the 
local inextensibility constraint r' 2 (s,t) = 1. Parametriz- 
ing the contour r{s,t) = {r±,s — r^) T by its transverse 
and longitudinal displacements from a straight line (see 
Fig. [T]), the weakly-bending limit of small contour gra- 
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dients r' 2 = 0(e) <C 1 is realized for very stiff poly- 
mers (e = L/£ p ), alternatively for semificxible filaments 

strongly prestretched with a force / prc (e = fp^J 2 /£ v ). 

The conformational dynamics in solution follows from 
a balance of elastic and tensile forces —SH/6r, thermal 
noise £, and anisotropic friction [r'r' + £(l — r'r')]<9 t r 
Within the weakly-bending limit, transverse and longi- 
tudinal fluctuations have strongly different correlation 
lengths: £x/£\\ = C(e 1/2 ); cf. Table Hfb). An adiabatic 
approximation (justified via a multiple scale analysis) ex- 
ploits this scale separation. The resulting equations of 
motion [l3| are written in terms of formally independent 
rapidly and slowly varying arclength parameters s and 
se 1 / 2 , respectively: 

d t r± = -d^rx + fd 2 rj_ 
d 2 s f = -C(d t g). 
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Eq. (|2"a|) gives the small-scale dynamics of the transverse 
displacements r±(s,t) for locally constant tension / = 
f(s, t), cf. Eq. |T]). Using a Cosine transform with respect 
to s, it is readily solved by the response function 

X x(q;t,t') = e -« 2 [« 2 (*- t ')+/ t 'dT/(,- ) r)] ^ _ t ,y (3) 



Eq. (|2b[) describes the coarse-grained tension variations 
on the large scale se 1 / 2 : it relates curvature in the tension 
to (average) changes in stored length density {g) (s, t) = 

_i_ / (s, t). Averaged both thermally and spatially (on 
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the small scale s), (g) inherits its remaining s-dcpcndence 
from the tension / in Eq. (|3|): 
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f°— f dt'q X ±(q;t,t')^(q,t') 
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(4) 



Reintroducing a single unique arclength variable, s = s, 
Eqs. (|2bj) and (j4|) result in a nonlinear partial integro- 
differential equation (P1DE) for f(s,t), that was ana- 
lyzed in Ref. [13J for explicitly prescribed boundary con- 
ditions. In the present case, however, the boundary con- 
dition at s = has to be determined implicitly. The 
polymer's incxtcnsibility requires that the bulge be cre- 
ated using stored length from the tails. To formalize 
this condition, we demand at any time a vanishing aver- 
age longitudinal velocity (9tf\\) at the origin where the 
force is applied, and also at infinity. Inextensibility (ri = 

\r'l + 0(e 2 ) « g) gives = f °°d S (9 t r|,) = J °°ds (d t g). 
With d s f\ s ^, 00 = and Eq. (|2bj) . this constraint implies 



d s f\ s =o = -( / dsd t (q - g) 



(5) 



The difference (g — ~g) represents the excess length stored 
in the bulge on the small length scale £± . Consequently, it 
did not contribute to Eq. which was spatially coarse- 
grained on intermediate scales £± -C I -C £\\. It can be 
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FIG. 2: (Color online) Numerical solutions f(s,t) to 
Eqs. (|2bl4l6p for / prc = 0, time is increasing from light to dark 
color. Inset: log-log plot of the effective longitudinal force 
/|| (t) = /(0) (circles/dashed), and of the nonlinear response 
Aj_(t) from Eq. (JJJ (solid). Dotted lines indicate the asymp- 
totes of Table|I|b). The crossover scales are ti = {jf±)~ 2 and 

»f = (^/C) 1/2 (7/±)- 1/4 , with 7 = (^C) 2/7 /l /7 - 



obtained, though, from the right hand side of Eq. Q 
upon replacing £± — > — /j_ sings 0(f). Evaluating the 
s-integral in Eq. ([5|) to leading order yields our central 
analytical result: a boundary condition for the tension 
that quantifies the feedback between "bulge" and "tail" 
dynamics: 



o F, C/I f°°dq 
d s f\ s =o = y- / —d t 

4 In 7T 



t -,2 

df qxx(q;t,t')\ s=0 



(6) 

In terms of the response function X x{q\ t, t') of Eq. j3|), 
the average displacement Aj_(£) induced by the trans- 
verse force (i.e., the nonlinear response) reads 



A±(t) = fx 



I — I dt'xx(q;t,t%= , 

Jo 71" Jo 



(I) 



which is evaluated at s = after the tension profiles 
f(s,t) are computed from Eqs. (12bl4l6l). To this end, 
we introduce two- variable scaling forms [13| that remove 
any parameter dependence: f(s,t) = jf±ip(s/sf,t/t{), 
with the crossover scales tf and Sf and 7 as in Fig. [5] 
Numerical solutions are obtained by mapping the PIDE 
onto a system of nonlinear equations [16[. Selected ten- 
sion profiles are displayed in Fig. [2] and describe one 
half of the filament with fj_ being applied at the ori- 
gin. Our analytical approach is based on reducing the 
scaling forms (p(s/s{,t/t{) to one-variable scaling func- 
tions ip ~ (i/if) Q ^(s/^||(£)) with £\\(t) = Sf{t/t { ) z in the 
asymptotic limits of short and long times. In the latter 
limit t 3> if, we recover either the taut-string approx- 
imation of Ref. [1| and may neglect bending and ther- 
mal forces, or the quasi-static approximation of Ref. [If3 |. 
which lets us treat the tension as locally equilibrated. 
Which approximation is valid depends quite strongly on 
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FIG. 3: (Color online) Regimes of intermediate asymp- 
totics (separated by thick black lines) for the nonlinear re- 
sponse Ax(t) (boxed formulas); time t/( , yf±)~ 2 vs. force ra- 
tio ,/pro/(7/±) (log- log scale). The universal initial regime 0] 
(light shaded) is followed by a quasi-static regime (white) with 
different force scaling for asymptotically small (<) and large 
(>) 01 force ratio; in these limits, the respective prefactors 
are b< ~ [8(1 + V2) 2 /tt s ] 1/s and b> ~ tt~ 1/2 . An intermediate 
taut-string regime (dark shaded) emerges for very small force 
ratio. The prefactor is a = [3(2 + \/2)/7r 2 ] 2/9 if / pro = 0. 



the prestretching force / pre through the ratio / pre / (7/1)7 
similar to the related scenario of longitudinal stretch- 
ing forces applied to prestretched filaments The 
resulting intermediate asymptotic scaling laws for A_l (t) 
are summarized in Fig.[3j including analytical prefactors. 
For a given ratio / pr c/(7/jJ, the evolution of Aj_(i) cor- 
responds to a vertical path through Fig. [3] The exact 
solutions quickly converge to these asymptotes, as shown 
in the inset of Fig. [2] for the limiting case / prc = 0. 

In summary, we argue that the coupling between trans- 
verse and longitudinal response affects not only single 
polymers, but also single crosslinks , crosslinked networks, 
and tensegrity structures [!, 0, E2, E3]- For complete- 
ness, we note that our self-consistent approach both for 
the heuristic "bulge" idea as well as for the systematic 
derivation of Eq. (|6|) applies only to the nonlinear [1 31 ] re- 
sponse on sufficiently small times t -C t^,t c . At v L , end 
effects become important, and at t c , the weakly-bending 
assumption breaks down: the contour gradients become 
large when Aj_ — (±- We find that t c > t L for ini- 
tially weakly-bending filaments (as those in the above 
discussed situations) Our analysis of the generic 



coupling mechanism is not constrained by the details of 
the relaxation regime t > (J (which is similar to the 
athermal case). 
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Since the weakly-bending assumption still holds at t = if , 
longitudinal friction is the only relevant nonlinearity and 



higher-order terms oc r± in Eq. pa[) are negligible 
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if f± > ip and /P^ ^ 



otherwise 



t c > t\. For f ± > £ a J 2 /L 7/2 (/ pr8 > £l/L 4 ), 4 falls into 
the taut-string (quasi-static) regime (cf. Fig. [3} . 



